Quintessence, super-quintessence and observable quantities in Brans-Dicke and 

non-minimally coupled theories 
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The different definitions for the equation of state of a non-minimally coupled scalar field that have 
been introduced in the literature are analyzed. Particular emphasis is made upon those features 
that could yield to an observable way of distinguishing non-minimally coupled theories from General 
Relativity, with the same or with alternate potentials. It is found that some earlier claims on that 
super-quintessence, a stage of super-accelerated expansion of the universe, is possible within realistic 
non-minimally coupled theories are the result of an arguable definition of the equation of state. In 
particular, it is shown that these previous results do not import any observable consequence, i.e. that 
the theories are observationally identical to General Relativity models and that super-quintessence 
is not more than a mathematical outcome. Finally, in the case of non-minimally coupled theories 
with coupling F = l + £<^> 2 and tracking potentials, it is shown that no super-quintessence is possible. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

During the last years, observations of distant Type la 
supernovae Q and CMB measurements M have shown 
that the universe is, most likely, undergoing a process of 
accelerated expansion. The widespread vision of the cos- 
mological model is, since then, a spatially flat low matter 
density universe. This implies that the total energy den- 
sity today is dominated by a contribution having nega- 
tive pressure (cosmological constant, or quintessence ||) 
which has just began to undertake the leading role in the 
right hand side of the Einstein field equations. 

The cosmological constant solution to this state of af- 
fairs appears not to be completely satisfactory (see for 
instance Q). Precise initial conditions should be given 
in order to solve the coincidence problem (why the vac- 
uum energy is dominating the energy density right now). 
Moreover, a fine-tuning problem appears, since a vac- 
uum energy density of order ~ 10~ 47 GeV 4 requires a new 
mass scale about 14 orders of magnitude smaller than the 
electroweak scale, having no a-priori reason to exist. In 
addition, the equation of state, pj p, for vacuum energy is 
exactly equal to —1, what at first sight appears as yet an- 
other value which is, in a dynamical setting, precisely set. 
Quintessence ||, and its derived models, being the main 
alternatives, are based on the existence of one or many 
scalar fields, which dynamically evolve together with all 
others components of the universe. The above-mentioned 
problems are alleviated within this framework. A sub- 
class of models, those having inverse power law poten- 
tials, present tracking solutions where a given amount of 
scalar field energy density can be reached starting from 
a large range of initial conditions (see for e.g. Refs. ||). 

The simplest models of quintessence are based on mini- 
mally coupled scalar fields. For a general potential V, the 
equation of state for those quintessence models is given 
by 



p _ f - 2V 
P 



2V' 



(1) 



and it can be easily proven that this expression is 
bounded to be within the range -1 < p/p < 1, unless 
of course one is willing to accept negative defined poten- 
tials. In the latter cases, the energy density becomes itself 
a negative quantity. For usual models of quintessence, 
then, it is clear that no super-acceleration can appear. 
The latter is a result of an extremely negative (< — 1) 
equation of state. This possible super-accelerated ex- 
pansion has been recently dubbed super-quintessence by 
several authors, e.g. although its consequences are 
being analyzed since some time before g. The mam rea- 
son supporting this interest is that current observational 
constraints are indeed compatible with, if not favoring, 
such values for the equation of state (see for instance Ref. 

0). 

Extended quintessence models are those in which the 
underlying theory of gravity contains a non-minimally 
coupled scalar field. It is this same scalar field which, 
apart from participating in the gravity sector of the the- 
ory, is enhanced by a potential to fulfill the role of normal 
quintessence. From a theoretically point of view, these 
ideas are appealing: it is the theory of gravity itself what 
provides the dynamically evolving, and currently dom- 
inating, field. Recent works on this area include those 
presented in Refs. §, |, § 0, ||, ||, (l| . We will have 
the opportunity to comment with much more detail on 
some of these works below. In addition, just to quote a 
few others in a so vastly covered topic, see the works of 
Rcf . . We would also like to remark that one of the 
first detailed analysis of a non-minimally coupled the- 
ory with a scalar field potential was made by Santos and 
Gregory ||l5| , years before the concept of quintessence 
was introduced. 

It has been claimed by many that a non-minimally 
coupled theory, like for instance Brans-Dicke gravity, can 
harbor super-quintessence solutions (e.g. || P). How- 
ever, there are different, and in most cases conflicting, 
definitions for the equation of state in these theories. 
Then, care should be exercised when analyzing the claims 
of the existence of super-quintessence solutions: in some 
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cases, they do not report either any physical import, be- 
cause the equation of state really is not more than a 
complex relationship between the field and its deriva- 
tive without any supporting conservation law, nor any 
observational consequence, because the amount of super- 
quintessence is so small that is far beyond any foreseen 
experiment. It is the aim of this paper to help in clari- 
fying these points, and to analyze, from an observational 
point of view, how non-minimally coupled theories differ- 
entiate from usual General Relativity in what concerns 
to quintessence and super-quintessence models. 

The rest of this work is presented as follows. In the 
following Section we comment on the energy conditions 
and the status of super-quintessence regarding them. 
Then, we introduce the gravity theories we are interested 
in. Section IV analyzes the case of Brans-Dicke gravity 
whereas Section V studies more general non-minimally 
coupled theories. A discussion and summary of the re- 
sults is given in Section VI. A brief Appendix discusses an 
alternative formulation of the theories of gravity, useful 
for numerical computations. 



II. THE ENERGY CONDITIONS 

For a Friedman- Robertson- Walker space-time and a di- 
agonal stress-energy tensor = (p, —p, —p, —p) with p 
being the energy density and p the pressure of the fluid, 
the energy conditions (EC) read: 



null: NEC 
weak: WEC 
strong: SEC 
dominant: DEC 



(p + P>0), 

(p>0) and (p + p> 0), 

(p + 3p > 0) and {p + p > 0), 

(p> 0) and (p±p> 0). (2) 



They are, then, linear relationships between the energy 
density and the pressure of the matter/fields generat- 
ing the space-time curvature. Violations of the EC have 
sometimes been presented as only being produced by un- 
physical stress-energy tensors. If NEC is violated, and 
then WEC is violated as well, negative energy densities 
-and so negative masses- are thus physically admitted. 
However, although the EC are widely used to prove the- 
orems concerning singularities and black hole thermody- 
namics, such as the area increase theorem, the topological 
censorship theorem, and the singularity theorem of stellar 
collapse fll6| ] , they lack a rigorous proof from fundamen- 
tal principles. Moreover, several situations in which they 
are violated are known; perhaps the most quoted being 

for 
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the Casimir effect, see for instance Rcfs. 
ditional discussion. Observed violations are produced by 
small quantum systems, resulting of the order of fi. It is 
currently far from clear whether there could be macro- 
scopic quantities of such an exotic, e.g. WEC-violating, 
matter/fields may exist in the universe. A program for 
imposing observational bounds (basically using gravita- 
tional micro and macrolensing) on the existence of matter 



violating some of the EC conditions has been already ini- 
tiated, and experiments are beginning to actively search 
for the predicted signatures [[18]. Wormhole solutions to 
the Einstein field equations, extensively studied in the 
last decade (see Refs. |ll| for particular examples), 
violate the energy conditions, particularly NEC. Worm- 
holes are probably the most interesting physical entity 
that could exist out of a macroscopic violation of the 
EC. 

It is interesting to analyze what does super- 
quintessence imply concerning the validity of the EC. 
As stated in the Introduction, super-quintessence is de- 
scribed by a cosmic equation of state 



?<-l, 
P 



(3) 



and so different situations arise depending on the sign of 
the energy density p. If p > 0, super-quintessence implies 
p + p < 0, and thus the violation of all the point-wise EC 
quoted above. Note that WEC is violated because of the 
violation of its second inequality. If, on the contrary, al- 
ready p < 0, then NEC may be sustained, but WEC is 
violated. Super-quintessence then implies strong viola- 
tions of the commonly cherished EC. But, should this be 
taken as sufficiently unphysical as to discard a priori the 
possibility of a super-accelerating phase of the universe? 

Apart from the finally relevant response, coming from 
experiments (today super-quintessence equations of state 
are not discarded, and maybe even favored by experimen- 
tal data, see for e.g. Q), the answer will of course rely 
on how much do we trust the EC, which, as we have al- 
ready said, are no more than conjectures. Particularly 
for non-minimally coupled theories, violations of the EC 
are much more common than in General Relativity, see 
for instance the works of Ref. and references therein. 
In addition, recently, the consequences of the energy con- 
ditions were confronted with possible values of the Hub- 
ble parameter and the gravitational redshifts of the old- 
est stars in the galactic halo ^0|. It was deduced that 
for the currently favored values of Hq, the strong en- 
ergy condition should have been violated sometime be- 
tween the formation of the oldest stars and the present 
epoch. SEC violation may or may not imply the viola- 
tion of the more basic EC, i.e. NEC and WEC, something 
that have been impossible yet to determine. In any case, 
super-quintessence could be a nice theoretical framework 
for explaining observational data opposing the EC. To 
the study of super-quintessence in non-minimally cou- 
pled theories, we devote the rest of this paper. 



III. GRAVITY THEORY 

In this section we shall present the general non- 
minimally coupled Lagrangian density given by 



-V(4>)+L f i uld ] 



(4) 
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Here, R is the Ricci scalar and units are chosen such 
that 8nG = 1. The functions lu(<P) and V((f>) specify the 
kinetic and potential scalar field energies, respectively. 
The Lagrangian Lfi ui d includes all the components but 
(j). The function / will be assumed to be of the form 



f((l>,R) = F(4>)R. 
Einstein equations from the general action 



are: 



2 1 ( w 
ti - — \p fluid + "2 <; 



V-3HF) , 



H = ~2F 



{p fluid + P fluid) + up + F - HF 



3Hcl>=-—{ uj^ z - F+R + 2V, 



(5) 



(0) 



(7) 



(8) 



where overdots denote normal time derivatives. The 
Klein-Gordon equation is actually very complicated in 
the general case. Using that 



R = 6(H+ 2H 2 



1 

F 



P fluid - Zp fluid - uj<]> 2 +4V- 3(F + 3HF) , (9) 



after some algebra, it ends up being 




= -ZH6 lo 

2lo ' 



-V, 



F, 



F 



{ — P fluid + 3p fluid) — 7^7 



loF 



LO 



2V 



2Flo 



3F, 
2uF' 



■F, 



F' 2 



2ujF 



9Hcf). 



(10) 



Two different kinds of theories are usually studied. 
One is the archetypical Brans-Dicke gravity |30|] , which 
appears by choosing lo = lobd/4> and F = (j), where 
lubd is referred to as the coupling parameter. The 
other, generically named as non-minimally coupled the- 
ories (although of course Brans-Dicke gravity also has a 
non-minimally coupled scalar held), are those for which 
lu = 1, and F and the potential V are generic functions 
of the held. Interesting differences appear when in the 
latter cases is of the form F = const. + g(4>), they will 
be discussed below. At least formally, starting from one 
of these Lagrangian densities, one can always rephrase it 
into the alternative form by a redefinition of the scalar 
field. Sometimes, however, this cannot be achieved with 
closed analytical formulae. 



that deviate from those of GR by less than this amount 
in the current cosmological era. In general, these devi- 
ations from GR can be specified by the post-Newtonian 
parameters p3 



(dF/d(/)) 2 



(3-1 



luF + {dF/d<j)) 2 



1 F{dF/d<t>) d-y 
\2ujF + Z{dF/d4>) 2 l4>' 



Solar system tests currently constrain p3| : 
|7-1|<2 1CT 3 , |/3 — 1| < 6 10 



(11) 



(12) 



(13) 



and they translate into a limit on 
the current time, supposing lo 



l/F{dF/dcj)) 2 at 
= 1, specifically 
\/F(dF/d<t>) 2 < 2 10- s If on the contrary, we as- 
sume the form of Brans-Dicke theory, they imply lubd > 
500. This value has been derived from timing experi- 
ments using the Viking space probe J2^| . In other situa- 
tions, claims have been made to increase this lower limit 
up to several thousands, see Ref. Q for a review. 

Starting from the action, one can define the cosmolog- 
ical gravitational constant as 1/F. This factor, however, 
does not have the same meaning than the Newton gravi- 
tational constant of GR. The Newtonian force measured 
in Cavendish-type experiments between two masses mi 
and rri2 separated a distance r is GeffWi^/'' 2 , where 
G off is given by §1 



1 {2uF + 4Fl 

(Jeff — -=: 



F \ 2ujF + 3Fj 



(14) 



The previous expression reduces to the well-known equal- 
ity G cS (t) = (l/0(i))(2u> + 4)/(2w + 3) for Brans-Dicke 
theory. Current constraints imply 



, G e ff 
I / > 

Lreff 



< 6 10- 12 yr _1 - 



(15) 



In general, though, one cannot make the statement that 
this constraint does directly translate into one for F/F, 
for one could in principle find a theory for which even 
when F varies significantly, G c g does not. Example of 
this is the case of Barker's theory [ p5[ , where G c g is 
strictly constant. 

Nucleosynthesis constraints can also be set for to, how- 
ever their impact is smaller than those set up in current 
experiments (see for e.g. Refs. |^(| and articles quoted 
therein). 



A. Experimental constraints 



B. The General Relativity limit 



The predictions of General Relativity in the weak field 
limit are confirmed within less than 1% f23fl . Any scalar- 
tensor gravity theory, then, should produce predictions 



There are important differences between Brans-Dicke 
gravity and more general non-minimally coupled theo- 
ries, particularly in what refers to quintessence. 
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In the case of Brans-Dicke, when the coupling param- 
eter is large, the field decouples from gravity, and the 
theory reduces itself to General Relativity jwj. When 
there is a potential, the limit of u>bd - > oo would make 
the theory GR + A for every V. This is certainly not the 
case in non-minimally coupled theories when F involves a 
term independent of the field (a constant). The limiting 
case of a non- variable -F-function is, in that situation, not 
GR plus a cosmological constant, but GR plus the same 
potential. In this case, then, the field recovers the status 
of normal quintessence, being it minimally coupled and 
enhanced by a generic potential. It is only in this sense 
that comparing different theories with the same potential 
is justified. The same procedure do not provide meaning- 
ful results when working with Brans-Dicke (or induced 
gravity) models. To see how this difference appears it 
would be enough to focus on the different Klein-Gordon 
equations for both theories. In the case of Brans-Dicke, 



(p-3p) 



' + 3-0 = 

a IijJbd 



Iujbd 



2V - 



dV 



(16) 



and all terms in the right hand side (including those hav- 
ing the potential) are proportional to X/u>bd- Then, a 
sufficiently large value of wbc will make this equation 
sourceless, i.e. a solution being = const., and reduce 
any V in the Lagragian to a constant as well. This does 
not happen in the case of non-minimally coupled theories 
where F contains an independent factor. For instance, in 
the case in which F = l+£0 2 , the Klein-Gordon equation 
is 



3H(f> = £i?0 



dV 



what clearly shows that the limit £ 
theory into normal quintessence. 



(17) 



converts the 



IV. BRANS-DICKE THEORY 
A. Field equations 

Consider the Brans-Dicke action given by 

d 4 x 



-9 



J fluid ■ 



(18) 

We shall consider that the matter content of the universe 
is composed by one or several (non-interacting) perfect 
fluids with stress energy tensor given by 



(19) 



where v^v^ = —1. This previous equation, then, with ad- 
equate values of p and p will be valid for the contributions 
of both, dust and radiation. Finally, we shall assume that 
the universe is isotropic, homogeneous, and spatially flat, 



and then represented by a k = Friedmann-Robertson- 
Walker model whose metric reads 



ds 2 = -dt 2 + a 2 (t)[dr 2 + r 2 d6 2 + r 2 sin 2 6d<j> 2 }. (20) 
In this setting, the field equations are given by 
- 2 d<j> uj4> 2 V 



a a <p to <p V p 
2 a 6 <j) 2 30 30 ' 



( a uj 2 V 
! a0 + 2 0^ ~ "0 



d- = (p-3 P ) ^_ 
a 2w + 3 2w + 3 



2V 



,dV 



(21) 



(22) 



(23) 



To simplify the notation in this Section we shall name 
ujbd simply to. The continuity equation follows from the 
Bianchi identity, yielding the usual relation 



p + 3-(p + p) =0, 



(24) 



which applied to both, matter (p m — 0) and radiation 
(p r = 1/3 p), gives the standard dependencies: 



Pm — Pm,Q a 



Pr — Pr,O a 



(25) 



We have chosen the scale factor normalization such that 
a at the present time is a = 1. The current values of the 
densities are given, in turn, by 



Pm,0 



3H 2 Q 



^ '771,0? 



Pr,o — 3H 2 il r fi. (26) 



Here, Ho = 100 x h km/s/Mpc is the current value of 
the Hubble parameter and fi r = 4.17 x 10 -5 //i 2 is the 
radiation contribution to the critical density (taking into 
account both photons and neutrinos, see pl| ). Typically, 
we shall work in a model with Q m = 0.4, but this can 
be fixed to any other value we wish, by using the contri- 
bution of the Brans-Dicke field to respect the flatness of 
the universe. 

The contribution of the field to the field equations can 
be directly read from the field equations, if we replace 
the usual General Relativity gravitational constant with 
the inverse of 0. The effective energy and pressure for 
the field end up being, 



V 
"3 



a 

— ( 

a 



(27) 



and 



- ^--V + + 2-c 

2ch a 



(28) 
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B. Numerical implementation 

After being unable to find any obvious coordi- 
nates/field transformation, in the sense explored by Mi- 
moso & Wands ^7j, Barrow & Mimoso and Barrow & 
Parsons (28), and Torres & Vucetich that can deal 
with the complexities introduced by the appearance of 
the self-interaction and solve the system analytically, we 
have prepared a computed code to integrate the system 
( pl[ - |2^ ) numerically. Indeed, not all 4 equations in the 
system are independent, because of the Bianchi identi- 
ties, and we have chosen to integrate Eqs. (|l|) and j2^ ) 
havingas input the form of the matter densities given in 
Eq. ( |25|) . We have followed the original idea of Brans- 
Dicke p(i|, [33[ , and transformed Eq. (pi| ) into an equation 
for H, by completing the binomial in the left hand side. 
Our variable of integration was x = ln(£), and the output 
were ln(et), (f> and <f>' , where a prime denotes derivative 
with respect to x. Having these values for each moment 
of the universe evolution, it is immediate to obtain </>, 
H , and any other quantity depending on them, like the 
effective pressure and density of the Brans-Dickc field 
given in Eqs. ( p7f]28| ). The relevant initial condition of 
the integration (in <fi) is chosen such that we fulfill to- 
day the observational constraint (G c g = 1) given by Eq. 
( |l4| ) . The derivative of the field can be set within a very 
large range at the beginning since, while producing unob- 
servable changes at the early stages of the universe, this 
initial condition is washed out by the evolution (a large 
range of different initial conditions will give the same re- 
sults). The potential is generically written as 

V = V = V f(<f>), (29) 

and the value of the constant Vq is iteratively chosen such 
that it fulfills the requirement of a large (say, fl^ ~ 0.6) 
field contribution to the critical density at the present 
time, for any given function f(<p). We have tested our 
code in the limiting cases of the problem and found agree- 
ment with previous results. As we have discussed, when 
lj — ► oo, Brans-Dicke theory becomes General Relativ- 
ity, (f> being a constant, and every potential effectively 
behave as a cosmological constant (i.e. v<t>l P4> = — 1) 
during all the universe evolution. Additionally, when we 
are in pure Brans-Dicke theory, without any potential, we 
reproduce the results of Mazumdar et al. for the ratio 
between the Hubble length at equality and the present 
one, a eq H eq /H 0. 

C. A worked example: V = Vo4>~ 2 

In Figure [j] we show the different contributions to the 
critical density during the universe evolution for a Brans- 
Dicke theory with u> = 500 and inverse square potential. 
The contribution of the field is given, at any time of the 
universe history, by 



T — I — I — I — I — I 
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FIG. 1: Evolution of the contributions to the critical density 
for an inverse square potential in a Brans-Dicke theory with 
to = 500. H = 75 km/s/Mpc and f2 m , = 0.4. 



and the others fi-values are defined in the same usual 
way as well. We see that the Brans-Dicke field can act as 
quintessence, in agreement with what other authors have 
previously found (see for example Refs. @, [l3| and ref- 
erences therein) . The value of Vq is extremely small, and 
mimics a cosmological constant in General Relativity. 
The Brans-Dicke field and its derivative do evolve in time. 
However, the current value of is <fi = 9.8 x 10~ 14 yr , 
fulfilling the above-mentioned constraint. The equiva- 
lence time (i.e. when p ma tter = Pradiation) in this model 
happens at 19801 yr, or ln(a) = —8.59. 

In Figure ||we show, for the same model, the evolution 
of the ratio between the effective pressure and densities 
for the Brans-Dicke field. This ratio evolves strongly 
during the recent matter era, the reason being that p$ 
actually crosses zero (from negative to positive values). 
This is in agreement with Figure [I], which shows the 
current field domination. 

Indeed, we can obtain similar results to those presented 
here changing the form of the potential to a variety of 
functional dependencies on the Brans-Dicke field (see Ta- 
ble I). Most interestingly, we see that at the present age 
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FIG. 2: Evolution of the ratio between the effective pressure 
and densities for the Brans-Dicke field (p^/p^), and of the 
effective equation of state for the field (w^>), entering in the 
observable quantity H . The latter quantity is defined in Sec- 
tion IID. Model parameters are as in Figure 1. 



of the universe, the effective equation of state for some 
potentials [we remark here that this is an abuse of lan- 
guage, as will be explained below] is smaller than -I. As 
we stated, the phenomenon of having p^/p^ < —1 has 
been referred to as 'super-quintessence' by other authors, 
whereas the corresponding dark energy has been dubbed 
'phantom energy' Apparently, already the simplest 
potentials one can imagine can be super-quintessence po- 
tentials within Brans-Dicke theory. It is true, however, 
that the amount of super-quintessence we have found 
(how large is deviation from -1 towards smaller values) 
is very small. It is indeed much smaller than what other 
authors claimed before (see in particular Ref. ||). How- 
ever, we note that, apparently, there is a sign mistake in 
their Klein-Gordon equation (7), the last term in their 
rhs should be positive, what can be tested by differenti- 
ation or comparison with, for example, Eq. (2.3) of Ref. 
jjJJ or Eq. (2.6) of J|2| This can actually produce a much 
bigger difference from an equation of state equal to -1 (as 
we numerically tested) , and is probably the origin of the 
discrepancy. 



As we have briefly implied above, p^/p^ < — 1 do not 
have a especially clear physical meaning. Both, and p^ 
are made up of terms coming from the Lagrangian density 
for the field. But they also contain terms coming from 
the interaction between the field and gravity (through its 
non-minimal coupling). The crucial aspect, then, is that 
the ratio p,p/p^, does not represent an equation of state, 
like those of the other components, since an equation of 
the form p^ + 3d/a(p c / ) + p$) = is not fulfilled. We can 
actually see from first principles why p^+^a/a{p^+p^) =^ 
. The field equations of Brans-Dicke theory in a general 
metric are 

G» v = ^{T™ atter +Tp, (31) 

where is the Einstein tensor, T™ atter is the stress- 
energy tensor for the matter sector of the theory and 

(32) 

with A being the D'Alambertian operator. Now, if we 
multiply the previous equation (|T]) by 6 and take covari- 
ant derivatives, it can be seen that |30| 

Tp * = m^). u (33) 

so that, whereas the usual continuity equation for matter 
is valid, the continuity equation for the above-defined 
'stress-energy tensor' for the field gets complicated. 

D. Observable quantities and the equation of state 

Then, if not p^/p^, which is the relevant (physically 
meaningful) quantity to be considered as the equation of 
state for the field in Brans-Dicke theory? We suggest that 
the important quantity to look at should come from what 
we actually measure. In the case of the homogeneous 
problem we are analyzing, this is the Hubble parameter, 
H . If we rewrite the first of the Friedmann-Robertson- 
Walker equations as 

H 2 = - (p fluid + P<t>) 

= I (p,n,oa- 3 + Prfi a- A + p^ a-^ 1+w ^) , (34) 

then, is what is going to establish the departure from 
the predictions of General Relativity plus cosmological 
constant or a generic quintessence potential of a mini- 
mally coupled field. To be specific, if w<p — —I, the the- 
ory would be indistinguishable from General Relativity 
plus cosmological constant (from an observational point 
of view). If w,), > — 1, then the theory would be indistin- 
guishable from a normal (minimally coupled) scalar field 
with a given potential. And finally, only if < — 1, 
the theory would be observationally different from its 
general relativistic counterparts: w$ < —1 is a value 
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TABLE I: Values of the effective and phenomenological 'equa- 
tions of state' for the Brans-Dicke field for different poten- 
tials. Model parameters are u = 500, Q. r = 4.17 X 10" 5 //i 2 , 
Q m = 0.4, and H = 75 km/s/Mpc. 
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\J?VI < rf / today 


\ wtouay 


flr 1/2 


-0.9967 


-0.9973 


1 / 2 


-0.9986 


-0.9993 




-0.9947 


-0.9953 


3 / 2 


-1.0006 


-1.0013 


0" 1 


-0.9957 


-0.9963 


2 


-1.0016 


-1.0023 


0- 2 


-0.9937 


-0.9943 


4 


-1.0056 


-1.0063 


0" 4 


-0.9897 


-0.9903 


6 


-1.0096 


-1.0103 


r 6 


-0.9857 


-0.9863 




-0.9997 


-1.0003 




-0.9957 


-0.9963 


cos(0) + 1 


-0.9966 


-0.9972 



that cannot be attained by any minimally coupled po- 
tential, as we discussed in the Introduction. In that case, 
super-quintessence adopts its rightful meaning: a super- 
accelerated expansion of the universe, contrary to the 
case in which it just represent a particular relationship 
between quantities importing no physically clear concept. 
In the previous equation, and to be consistent with the 
generalized field equations of the theory, we have defined 



P — P0 a — —Pm,O a 



PrflO. 



{PmfiO, 



-3 



Pr,oa 



(35) 



where sub-indices m, r and stand for matter, radiation, 
and current values, respectively. By definition, 



3 \P<t>oJ m (a) 



(36) 



In Figure we have already shown the evolution of this 
quantity in time (redshift). We also see that, although 
some non-minimally coupled Brans-Dicke field can ef- 
fectively produce super-quintessence when supported by 
particular potentials, the effective 'equation of state' be- 
ing less than -1, its value is too close to -1 so as to 
effectively mimics a cosmological constant at the cur- 
rent epoch. In Table I we show that the same, al- 
most imperceptible, deviation from an effective cosmo- 
logical constant appears for other potentials. This is 
then showing a sign of caution when analyzing the impact 
of non-minimally coupled theories using the equation of 
state: differences are actually too small and fall below 
the threshold of any current or foreseen experiment. 



E. The meaning of 

Is Wtj, representing an equation of state in the usual 
sense? In other words, is the equation p^ — —3H(1 + 
w^pcj, , satisfied? To answer this question we shall 
rewrite the field equations for a general non-minimally 
coupled theory as [fL0| 



G 



flis 



(37) 



where we have defined a new stress-energy momentum 
tensor, on purpose, to make the previous equation valid: 



rp rjifluid 



(38) 



As a consequence of the contracted Bianchi identities T^ v 
is conserved. And since there are no explicit coupling 
between matter and fields, their corresponding energy- 
momentum tensors are also separately conserved: 



WT fluid = V»%u[<t>] = 



(39) 



In this framework, reporting no more than a rewriting of 
the field equations i.e. introducing no new physics at all, 
the explicit expression for T M „[0] is given by 



Vg^u + V M V„F - g^AF + (1 - F) G M „ 



(40) 



In order to obtain Eq. (j37|) we have just added and sub- 
tracted G^ v on the left hand side of the generalized Ein- 
stein equations, i.e. 



Grp rpfluid , rpfield 

with the usually defined 

1 



rpfield 
1 p.v 



V M V^ - g^AF (41) 



This immediately fixes the expression of T^ v [<f)\ above. 
The new effective energy density and pressure that this 
defined stress-energy density produces are 



ptp = lu^y + V(<J)) — 3HF + 3H (1 — F) , (42) 



p tk = L0-^--V((p)+2HF+F-(3H J +2H)(l-F) , (43) 

In the case of Brans-Dicke theory, we recall, u> — 
lu bd /(/} and F = 4>. Note then that p^ = p^ + 3H 2 {l — F), 
where p^ was given in Eq. (^7|). The defined equa- 
tion of state, thus, is exactly that given by p^/p^,, 
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since it was defined using the same Friedmann equation 
H 2 = l/3(pfiui d + Pd>)- Indeed, 

P<jj — —P fluid + ~p(p fluid + P<}>) 



— —P fluid + 3H 2 
= -(3H 2 F- p^) + 3H 2 
F) 



3H 2 (1 



At the same time, we can see that 



(3H 2 + 2H) (1 - F) . 



(44) 



(45) 



We conclude that the definition for represents a real 
equation of state, since it is supported by a conservation 
law, and that it is this the one that should be taken into 
account to compare with the predictions of GR, since it 
is directly related to the observable, H . We can also see, 
from Table I, that the difference between and the 
ratio P4>/ P(j> is very small. The reasons that leads to this 
are explicitly discussed for NMC theories in Section III, 
a similar argument applies here as well. 



F. CMB-related observables 

The evolution of the comoving distance from the sur- 
face of last scattering (defined as z = 1000, equivalently 
In a = —6.90) can be computed, for different theories, as: 



dr = 



da 



a 2 H{a)' 



(46) 



Only in the case of extremely low coupling factors (e.g. u) 
of Brans-Dicke theory), discarded by current constraints, 
we see a noticeable difference with the result of General 
Relativity plus cosmological constant. To give a quanti- 
tative idea we can quote the ratio 

tbd - tor 

T GR 

calculated today (a — 1), which, for to = 25 results equal 
to -0.014, whereas for u> = 500 is —5.7 10" 4 , and quickly 
tends to zero for bigger values of w. 

The angular scales at which acoustic oscillations occur 
are directly proportional to the size of the CMB sound 
horizon at decoupling, that in comoving coordinates is 
roughly and inversely proportional to the co- 

moving distance covered by CMB photons from last scat- 
tering until observation, that is tq — Td ec |34f| . The rmil- 
tipoles scale as the inverse of the corresponding angular 
scale, and so 



TQ — T de 
7~dec 



(47) 



As in the non-minimally coupled models studied in Ref. 
p6| , t changes because of a different dependence of the 



Hubble length H^ 1 in the past. However, we have al- 
ready noted that this change is almost imperceptible 
when compared with usual General Relativity plus a cos- 
mological constant, unless of course (violating current 
constraints) the coupling parameter uj is low enough. 

The Integrated Sachs- Wolfe effect makes the CMB co- 
efficients on large scales, small £'s, change with the varia- 
tion of the gravitational potential along the CMB photon 
trajectories [p9| . This is undoubtedly changed because of 
a variation in the gravitational constant since the time 
of decoupling. However, we expect this change to be also 
very small, since the G-variation we have found, for val- 
ues of uj = 500 and bigger, are typically less than 2% 
since the time of decoupling. 

The scale entering the Hubble horizon at the matter- 
radiation equivalence is also important, since it will de- 
fine the matter power turnover p5fl . The shift in the 
power spectrum turnover is given by |3q] 



5k t 



f SH- 1 



(48) 



and again, this reports a very small difference for all cur- 
rently possible w-values. Only for to = 25 this difference 
is about 12% (where a case of a power law potential with 
exponent equal to -2 taking as an example). For u> = 500 
and bigger, the differences are less than 1% (to give a 
precise example it reports a 0.7% difference in the same 
power law case as commented before and uj — 500). Con- 
trary to what Baccigalupi et al. have done in the past 
H , we are not comparing two different theories (Brans- 
Dicke and General Relativity) with the same potential, 
but rather, and motivated by the findings of this section, 
Brans-Dicke theory with any given potential and General 
Relativity plus A. It is in this case that the possibilities 
of actually distinguishing both situations diminish. 



V. NON-MINIMALLY COUPLED THEORIES 

The details of the cosmological evolution using the gen- 
eral action (0) above, with lu = 1 and 



F(4>) = i + t4> 2 



(49) 



were explored in Refs. || |T(j, [36|, among many others, 
and we refer the reader to these works and references 
therein for additional relevant discussions. Our nu- 
merical code is in agreement with the results therein 
reported, and it is a direct extension of the numerical 
implementation reported in Section IIIB. 

In this Section, following the previous discussion, we 
would like to focus on the possible definitions of equations 
of state and their impact onto observable quantities. Just 
as an example, we show in Figure |^ the case of a tracking 
potential of the form V = Vq<\>~ 2 , Hq = 70km/s/Mpc, in 
a flat universe with f2 m ,o = 0.4. The equivalent Brans- 
Dicke parameter (obtained by redefining fields in Eq. 
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FIG. 3: Matter, radiation and quintessence energy densities, 
both non- minimally (NMC) and minimally coupled (MC). 
The incomplete curve corresponds to the non-minimally cou- 
pled quintessence, and the cutoff is produced by a change in 
the sign of the effective density, see text for details. 



(||) in order for it to look like a Brans-Dicke theory) is 
lojbd = Y F ^, and its value is given by defining £. The 
value of £ used in the model of Figure || and successive 
ones is 5.8 10 -3 , what implies an equivalent Brans-Dicke 
parameter lobd = 3071, well in agreement with current 
constraints. 

Starting from the general field equations, we can 
immediately define, as done for the Brans-Dicke the- 
ory, an effective energy density and pressure for the 
scalar field. The former, for instance, appears writ- 
ing the 00-component of the field equation as H 2 = 
1/3F (p fi U id + ptj)) . The explicit expression then being, 



P4> = 2^ 

for the energy density, and 

P<t> = 2<P - 



V-3HF 



V + 2HF + F 



(50) 



(51) 



for the pressure. These two expressions do not, as 
we have shown before, pertain to a conserved energy- 
momentum tensor. They do, however, define an effective 



equation of state, this being just p^/p^. This relation- 
ship is subject to same caveats mentioned above for the 
case of Brans-Dicke: it is neither positive nor negative 
defined, since the effective energy density itself shifts its 
sign during the evolution. The energy density quoted 
above is what is depicted (whenever possible) in Figure 
|. As it can be seen, it tracks closely the usually defined 
minimally coupled (MC) energy density at low redshifts, 
this being an effect of the necessarily small coupling £ 
that is adopted to fulfill observational constraints. 

Again, in order to work with a conserved energy- 
momentum tensor, we can rewrite the field equations and 
obtain a real equation of state, = p^/p^, where p^ and 
P0 were given in Eqs. (ji|) and ([43]), respectively. As we 
already mentioned in the case of Brans-Dicke, this equa- 
tion of state is exactly what results in writing the field 
equation as H 2 = 1/3 (p fluid + Q<f>), defining implicitly 
P0 = 1/F(pf hu d + p<t>) - P fluid- Finally, just for compar- 
ison, one can as well consider the equation of state for 
the case of a minimally coupled scalar field, pj p, with 



p = \i> 2 + V(<t>) , 



p=\i> 2 - V (<P) 



(52) 



(53) 



In Figure [| we show the results of these different def- 
initions of equation of state for the current time, as 
a function of the exponent of the power law potential 
V = Vo <j) a , for a flat universe model given by fl m ,o = 0.4, 
the same model used in Figure ||. We see that they do 
not present noticeable differences. Very low values of 
the power law exponent (shallow potentials) are needed 
to produce equations of state near that generated by a 
cosmological constant. To be quantitative, Figure ^ also 
shows, in the right panel, the differences between the 
equation of state directly related with the Hubble pa- 
rameter, Wrf, and the non-minimally coupled (NMC) and 
minimally coupled p^/p^. Clearly, the differences are mi- 
nor. One can actually understand why these differences 
are so small. Note that the energy density and pressure 
in a non-minimally coupled field theory can be written as 
the minimally coupled ones plus additional terms. In the 
case of p<f,, these terms are equal to —3HF + 3H 2 (l — F), 
whereas for p^ only the first term above enters. Both 
these terms are, however, proportional to £, being them- 
selves 

-ZHF + 3H 2 (l -F) 6£H<j) - 3£H 2 cj) 2 . 

But since from the evolution of the field, 0(i7 _1 (/>)=f to- 
day, and at the current era, 0(H 2 ) ~ V, the energy den- 
sity can be written as p^ ~ p mc — 3£(^j-+F</> 2 ), where p mc 
is the minimally coupled energy density. Clearly, at the 
current cosmological era and because of the constraints 
on £, the second terms are sub-dominant in comparison 
to the first ones. A similar analysis can be established 
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FIG. 4: Left: Current values of the defined 'equations of state', see text, as a function of the power law exponent of the 
tracking potential. Right: Difference between the equation of state directly related with the Hubble parameter, w$ and the 
non-minimally coupled (NMC) and minimally coupled p^/p^,. 



for the pressure, where again all extra terms are propor- 
tional to £, and then for the equations of state. Today 
the influence of the last terms in Eqs. (^) and (|4 3| ), the 
'gravitational dragging' terms, as dubbed in Ref. |(l0||, is 
negligible in comparison to the minimally coupled con- 
tribution. It is only in the past, when the matter density 
dominates the evolution of the universe, that these terms 
become important. 

The previous analysis is not valid in the past history of 
the universe. Figure || shows the evolution of the different 
equations of state with redshift. We can see that the scal- 
ing solution of the tracking potentials (equation of state 
is approximately equal to — 2/(2+a) between z ~ 1-1000, 
when pmatter > P<j>, see for e.g. Refs. p| [To|), appears for 
all defined equations of state but w$. We also note that 
because of the sign non-definiteness of the defined non- 
minimally coupled effective 'energy density', the corre- 
sponding equation of state shows sudden changes at high 
redshift, in the position where the density cross zero. Fi- 
nally, we find that for tracking non-minimally coupled 
quintessence, there is no super-quintessence potential in 
any of the defined equations of state, not even by an small 
amount, all of them being greater than -1 at present. This 



latter result is actually valid for all potentials analyzed 
(similar cases than those presented in Tabic I). 

A. Comments on perturbations and on the possible 
degeneracy with kinetically driven quintessence 

Perturbations analysis of the held equations, as done 
by Perrotta and Baccigalupi JTcJ, have shown that rich 
phenomena are uncovered when working with the held 
equations written as in (|37|). The most important of them 
are those generated by gravitational dragging, exclusive 
of non-minimally coupled theories. This phenomenon is, 
basically, the early dominance of the last term in Eqs. 
( [42| ) and (f42|), which is, in turn, produced due to the fact 
that they are proportional to the square of the Hubble 
parameter, and then to the total energy density. The 
latter scales as a~ 3 and a~ 4 at early times, and then 
makes the gravitational generated term to dominate the 
dynamics of the held. 

The advantage of writing the formalism as in Eq. ( |3^ ) 
is the fact that usual perturbation analysis -as applied 
for any fluid component of the universe- is also valid for 
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be much different for non-minimally coupled gravity. In 
that case, 



o5 



in 
o 



cr 

CD 

"D 
CD 

H 

CD 
Q 




c eff . 



(55) 



for values of k ^> H . But because of gravitational drag- 
ging, 5p<f> can be quite difference from the usual case, and 
this ratio may be much lower than unity whenever the 
energy density perturbations of the scalar field are en- 
hanced by perturbations in the matter field. At the level 
of perturbations, then, quite distinctive effects appears 
in non-minimally coupled quintessence as compared with 
the usual case and make these theories possibly distin- 
guishable. 

Very recently, yet another scenario for an alternative 
model of quintessence was introduced p7] . In it, known 
as fc-essence, the Lagrangian density includes a non- 
canonic kinetic term: 



S = 



rf 4 :, 



1 



16nG 



R + P {4>^4>) 



Sn 



(56) 



where S m ^ r denotes the action for matter and/or ra- 
diation. Examples in which the Lagrangian depends 
only on the scalar field and its derivative squared 
X = -iV^V^. have been constructed |37j. The field 
equations for this models are 



FIG. 5: Evolution in redshift of the equations of state for a 
model with V oc <b~ 2 . 



(57) 



the field. In that sense, concepts as the equation of state, 
or the adiabatic sound speed 



1 



-<t> 



— = w <t> 



3H 1 



(54) 



(where we have made use of the field conservation equa- 
tion) can be well defined. If the field is slowly varying 
in time, ~ w$. Hu |3^] showed, for negative equa- 
tions of state, that adiabatic fluctuations are unable to 
give support against gravitational collapse. Density per- 
turbations would become non-linear after entering the 
horizon, unless the entropic term w^T^ > 39 1. 



The effective sound speed, cjg ^, is then defined in the 
rest frame of the scalar field, where ST^i = . The 
gauge invariant entropic term is written as w^T^ = (c^ ff — 
c %)^ eSt ^ > where 8^ est "> is the density contrast in the dark 

energy rest frame |nj S^ est) = (1+ w<f,) (v<j> - B) 

(we refer the reader to Ref. jlO| for more careful expla- 
nations). For normal quintessence, the effective sound 
speed is c^ s * = 1, giving a relativistic behavior to the 
corresponding density fluctuations. However, Perrotta 
and Baccigalupi [fLOl have found that the situation can 



where T™ ,r is the energy-momentum tensor for usual 
matter fields. p(<fr,X) corresponds to the pressure of 
the scalar field, whereas the energy density is given by 
Ptf, = 2Xdp/dX — p p8j . It can be seen that for this 
models, the speed of sound is given by M 



P,x 



-off 



PQ.X ^ 

PQ,X p,x + 2Xp : xx 



(58) 



Apparently, then, and since p is completely generic, it 
could exist a non-canonic kinetic term within fc-essence 
giving rise to the same results of non-minimally coupled 
gravity. Viscosity (a parameter relating velocity and met- 
ric shear with anisotropic stress |39| ) can however pro- 
vide the way to break the degeneracy, since it results 
non-zero for a non-minimally coupled field (contrary as 
well to what results in usual quintessence) |l0|. 



VI. CONCLUSIONS 

In the case of Brans-Dicke theory, and in the cases 
of WfiB allowed by current constraints, we have numer- 
ically proven that the homogeneous field equations of 
extended quintessence yield to no observable effect that 
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can distinguish the theory from the predictions of Gen- 
eral Relativity plus a cosmological constant. It is with 
this model that the comparison should be made, since 
for the large values of the coupling parameters required 
by current experiments, all potentials are closely similar 
to a constant function, and the theory itself to General 
Relativity. Although we have not made a detailed per- 
turbation analysis using the full numerical CMBFAST 
code, we can safely predict that the same situation will 
happen there, as a result of the analysis made for the 
CMB-related observables in Section IIC. We discussed 
the observationally-related definition of equation of state 
Wtf,, and not to the usually studied ratio between the ef- 
fective pressure and density directly obtained from the 
field equations, to which we assign no particular physi- 
cal meaning. The phenomenon of super-quintessence, i.e. 
a super-accelerated expansion of the universe, although 
possible for a non-minimally coupled Brans-Dicke scalar 
field, and impossible in any minimally coupled field situ- 
ation, it is of such an small amount that is far beyond the 
expectations of any realistic experiment. From a practi- 
cal point of view, then, it will always exist a scalar field 
potential supporting a minimally coupled field that pro- 
duces experimentally indistinguishable results from those 
obtained within the extended quintessence framework of 
Brans-Dicke theory. 

For the more general extended non-minimally cou- 
pled models studied, the possibility of having super- 
quintessence actually disappears: all tracking potentials 
explored produce effective equations of state greater than 
-1. We have shown that for low values of the exponent 
in the tracking potentials supporting the non-minimally 
coupled field (i.e. equations of state are close to -1), the 
difference among all defined equations of state is negli- 
gible. It is, however, in the perturbation regime where 
differences with the usual quintessence case can be no- 
ticed. As Perrotta and Baccigalupi have found juj, a new 
gravitational dragging effect appears here, giving rise to 
the possibility of clumps of scalar field matter. In this 
case, however, it is with fc-essence models that a degen- 
eracy could appear, particularly in those cases in which 
p.x + 2Xp^xx ^> 1, yielding the speed of sound to values 
close to zero. 

Finally, we remark that expanding solutions where 
acceleration is transient have been recently considered 
given the consistency problem between string theory and 
spaces with future horizons E3| . Since scalar-tensor theo- 
ries of gravity likely originate in string theory, it would be 
interesting to make a similar analysis to that presented 
in the previous reference for the case of non-minimally 
coupled theories. 
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Appendix 

In the literature, one may find two alternative in- 
troductions of general non-minimally coupled theories. 
Firstly, the one that we follow in Section II (see for in- 
stance p2|), and secondly, the one that is derived from 
the action 

i-/(^,i?)-^V^V^ 

-V(cf>) + L fluid \ , (59) 

where k is a constant, not necessarily taken as 1, and 
plays the role of the "bare" gravitational constant (see 
for instance jL2|, by the same authors). The function / 
is then assumed to be of the form 



{l/K)f{4>,R)=F((j>) R. 



(60) 



The function F, in the case we are interested in, is written 

as 



F(</>) = - + £0 2 

K 



i + e(0 2 -^)- 



(61) 



Then, a value of k = 1, as we have taken in the theoret- 
ical development of the previous sections just reduces to 
take 4>o, the current value of the field, equal to 0. This, 
however, is not what may result, in general, numerically 
convenient, since it would imply to precisely fix the evo- 
lution of the field to reach <fi = today. Instead, as we 
are not actually interested in any value of k per se, we do 
not make K = 1 in our numerical code, and instead follow 
the treatment given by Perrotta and Baccigalupi (|Tof ). 
In that case, they choose to rewrite the field equations as 



(62) 



with the corresponding field energy density and pressure 
given by 

P0 = / + V{(j>) - 3HF + 3H 2 (~ - f) , (63) 



wy -V{4>)+2HF - 



F - (3H 2 



2H) F 

1 K 



(64) 

When comparing Eq. (62) with the usual Einstein equa- 
tions, one has to take into account that the value of n is 
not 1 (although certainly it is truly close to unity, because 
of the constraint imposed on £). Then, if we decide to 



write the Friedman equation like H = l/3(pji U id + 
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to directly compare with GR (and the same matter den- 
sity) the corresponding relationship between and 
give in Eq. is 



(k - l)pfluid + KP4>- 



(65) 



In this scheme, p<f,/p<f, (with quantities defined as in Eqs. 
( |63| ) and (|64|)) will differ from the equation 



because Q4, ^ Pcj,, but it is the latter Eq. (66) what should 
be used to compare with the results of General Relativity 
with a fixed current matter density. 



1 1 1 ( e<t> 
o \6<po 



1 



\n(a) ' 



(66) 
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